Astro arXiv | astro-ph.HE

Does the gamma-ray binary LS I +61°303 harbor a magnetar? by Arthur G. Suvorov et al. on Wednesday 30 November The high-mass X-ray binary LS I +61{\deg}303 is also cataloged as a gamma-ray binary as a result of frequent outbursts at TeV photon energies. The system has released two soft-gamma flares in the past, suggesting a magnetar interpretation for the compact primary. This inference has recently gained significant traction following the discovery of transient radio pulses, detected in some orbital phases from the system, as the measured rotation and tentative spin-down rates imply a polar magnetic field strength of $B_p \gtrsim 10^{14}\,\mbox{G}$ if the star is decelerating via magnetic dipole braking. In this paper, we scrutinize magnetic field estimates for the primary in LS I +61{\deg}303 by analyzing the compatibility of available data with the system's accretion dynamics, spin evolution, age limits, gamma-ray emissions, and radio pulsar activation. We find that the neutron star's age and spin evolution are theoretically difficult to reconcile unless a strong propeller torque is in operation. This torque could be responsible for the bulk of even the maximum allowed spin-down, potentially weakening the inferred magnetic field by more than an order of magnitude. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2210.09471v2

Ultra high energy cosmic rays from past activity of Andromeda galaxy by V. N. Zirakashvili et al. on Wednesday 30 November It is shown that the relativistic jets associated with the growth and past activity of the supermassive black hole in the Andromeda galaxy could be the main source of cosmic rays with energies above $10^{15}$ eV. Most of the cosmic ray energy is related to a bow shock of the jet that produces multi-PeV cosmic rays with light composition. The highest energy cosmic rays with heavy composition are produced in the jet itself. The spectra of energetic particles produced in Andromeda galaxy and propagated to the Earth are calculated and compared with observations. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.04522v2

MAGIC observations provide compelling evidence of the hadronic multi-TeV emission from the putative PeVatron SNR G106 3+2 7 by MAGIC Collaboration et al. on Wednesday 30 November The SNR G106.3+2.7, detected at 1--100 TeV energies by different $\gamma$-ray facilities, is one of the most promising PeVatron candidates. This SNR has a cometary shape which can be divided into a head and a tail region with different physical conditions. However, it is not identified in which region the 100 TeV emission is produced due to the limited position accuracy and/or angular resolution of existing observational data. Additionally, it remains unclear whether the origin of the $\gamma$-ray emission is leptonic or hadronic. With the better angular resolution provided by these new MAGIC data compared to earlier $\gamma$-ray datasets, we aim to reveal the acceleration site of PeV particles and the emission mechanism by resolving the SNR G106.3+2.7 with 0.1$^\circ$ resolution at TeV energies. We detected extended $\gamma$-ray emission spatially coincident with the radio continuum emission at the head and tail of SNR G106.3+2.7. The fact that we detected a significant $\gamma$-ray emission with energies above 6.0 TeV from the tail region only suggests that the emissions above 10 TeV, detected with air shower experiments (Milagro, HAWC, Tibet AS$\gamma$ and LHAASO), are emitted only from the SNR tail. Under this assumption, the multi-wavelength spectrum of the head region can be explained with either hadronic or leptonic models, while the leptonic model for the tail region is in contradiction with the emission above 10 TeV and X-rays. In contrast, the hadronic model could reproduce the observed spectrum at the tail by assuming a proton spectrum with a cutoff energy of $\sim 1$ PeV for the tail region. Such a high energy emission in this middle-aged SNR (4--10 kyr) can be explained by considering the scenario that protons escaping from the SNR in the past interact with surrounding dense gases at present. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15321v1

Status and prospects of the CORSIKA 8 air shower simulation framework by Alexander Sandrock. on Wednesday 30 November The Fortran-versions of the CORSIKA air shower simulation code have been at the core of simulations for many astroparticle physics experiments for the last 30 years. Having grown over decades into an ever more complex software, maintainability of CORSIKA has become increasingly difficult, though its performance is still excellent. In 2018, therefore a complete rewrite of CORSIKA has begun in modern modular C++. Today, CORSIKA 8 has reached important milestones with a full-fledged implementation of both the hadronic and electromagnetic cascades, the ability to simulate radio and Cherenkov-light emission from air showers and an unprecedented flexibility to configure simulation media and their geometries. This presentation will discuss the current status of CORSIKA 8, highlight the new possibilities already available, and future prospects of this new air shower simulation framework. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16057v1

Continuous gravitational wave emission from neutron stars with pinned superfluids in the core by Brynmor Haskell et al. on Wednesday 30 November We investigate the effect of a pinned superfluid component on the gravitational wave emission of a rotating neutron star. Pinning of superfluid vortices to the flux-tubes in the outer core (where the protons are likely to form a type-II superconductor) is a possible mechanism to sustain long-lived and non-axisymmetric neutron currents in the interior, that break the axial symmetry of the unperturbed hydrostatic configuration. We consider pinning-induced perturbations to a stationary corotating configuration, and determine upper limits on the strength of gravitational wave emission due to the pinning of vortices with a strong toroidal magnetic field of the kind predicted by recent magneto-hydrodynamic simulations of neutron star interiors. We estimate the contributions to gravitational wave emission from both the mass and current multipole generated by the pinned vorticity in the outer core, and find that the mass quadrupole can be large enough for gravitational waves to provide the dominant spindown torque in millisecond pulsars. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15507v1

A Bayesian approach for torque modelling of BeXRB pulsars with application to super-Eddington accretors by A. S. Karaferias et al. on Wednesday 30 November In this study we present a method to estimate posterior distributions for standard accretion torque model parameters and binary orbital parameters for X-ray binaries using a nested sampling algorithm for Bayesian Parameter Estimation. We study the spin evolution of two Be X-ray binary systems in the Magellanic Clouds, RX J0520.5-6932 and RX J0209-7427, during major outbursts, in which they surpassed the Eddington-limit. Moreover, we apply our method to the recently discovered Swift J0243.6+6124; the only known Galactic pulsating ultra-luminous X-ray source. This is an excellent candidate for studying the disc evolution at super-Eddington accretion rates, for its luminosity span several orders of magnitude during its outburst, with a maximum $L_{\rm X}$ that exceeded the Eddington limit by a factor of $\sim 10$. Our method, when applied to RX J0520.5-6932 and RX J0209-7427, is able to identify the more favourable torque model for each system, while yielding meaningful ranges for the NS and orbital parameters. Our analysis for Swift J0243.6+6124 illustrates that, contrary to the standard torque model predictions, the magnetospheric radius and the Alfv\'en radius are not proportional to each other when surpassing the Eddington limit. Reported distance estimates of this source range between 5 and 7 kpc. Smaller distances require non-typical neutron star properties (i.e. mass and radius) and possibly lower radiative efficiency of the accretion column. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16079v1

Diffuse Emission of Galactic High-Energy Neutrinos from a Global Fit of Cosmic Rays by Georg Schwefer et al. on Wednesday 30 November In the standard picture of galactic cosmic rays, a diffuse flux of high-energy gamma-rays and neutrinos is produced from inelastic collisions of cosmic ray nuclei with the interstellar gas. The neutrino flux is a guaranteed signal for high-energy neutrino observatories such as IceCube, but has not been found yet. Experimental searches for this flux constitute an important test of the standard picture of galactic cosmic rays. Both the observation and non-observation would allow important implications for the physics of cosmic ray acceleration and transport. We present DINECRAFT, a new model of galactic diffuse high-energy gamma-rays and neutrinos, fitted to recent cosmic ray data from AMS-02, DAMPE, IceTop as well as KASCADE. We quantify the uncertainties for the predicted emission from the cosmic ray model, but also from the choice of source distribution, gas maps and cross-sections. We consider the possibility of a contribution from unresolved sources. Our model predictions exhibit significant deviations from older models. Our fiducial model is available at https://doi.org/10.5281/zenodo.7373010 . arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15607v1

MOJAVE XX Decade-long linear polarization variability in AGN jets at parsec scales by D. I. Zobnina et al. on Wednesday 30 November We studied the variability properties of the linear polarization of active galactic nucleus (AGN) jets on parsec-scales using stacking. Our sample is drawn from the MOJAVE program and consists of 436 AGNs manifesting core-jet morphology and having at least five VLBA observing epochs at 15 GHz from January 1996 through August 2019, with some additional archival VLBA data reduced by us. We employed a stacking procedure and constructed maps of (i) standard deviation of fractional polarization and electric vector position angle (EVPA) over epochs as the measure of variability, (ii) median polarization degree to quantify typical values in time. The distributions of these values were analysed along and across the jet. We found that EVPA variability in the core is typically higher than in the jet, most likely due to changes in opacity and component blending in the core region. The EVPA becomes more stable down the outflow. Most of the sources showing this trend have a time coverage of more than 12 years and at least 15 epochs. The possible cause could be the increase of stability in the magnetic field direction. The majority of AGNs exhibit insignificant trends of the relative fractional polarization variability along the ridgeline or across the jet width. There are no significant optical-class-dependent or spectral-class-dependent relations in the properties of EVPA and relative fractional polarization variability. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15624v1

The Rapidly Spinning Intermediate-Mass Black Hole 3XMM J150052 0+015452 by Z. Cao et al. on Wednesday 30 November A star tidally disrupted by a black hole can form an accretion disc with a super-Eddington mass accretion rate; the X-ray emission produced by the inner disc provides constraints on the black hole mass $M_\bullet$ and dimensionless spin parameter $a_\bullet$. Previous studies have suggested that the $M_\bullet$ responsible for the tidal disruption event 3XMM J150052.0+015452 (hereafter J150052) is $\sim$10$^{5} M_{\odot}$, in the intermediate black hole (IMBH) regime. Fitting multi-epoch XMM-Newton and Chandra X-ray spectra obtained after 2008 during the source's decade-long decay, with our latest slim accretion disc model gives $M_\bullet = 2.0^{+1.0}_{-0.3}\times10^{5} M_{\odot}$ (at 68% confidence) and $a_\bullet > 0.97$ (a 84.1% confidence lower limit). The spectra obtained between 2008-2014 are significantly harder than those after 2014, an evolution that can be well explained by including the effects of inverse-Comptonisation by a corona on the early-time spectra. The corona is present when the source accretion rate is super-Eddington, while there is no evidence for its effect in data obtained after 2014, when the mass accretion rate is around the Eddington-limit. Based on our spectral study, we infer that the corona is optically thick and warm ($kT_e=2.3^{+2.7}_{-0.8}$ keV). Our mass and spin measurements of J150052 confirm it as an IMBH and point to a rapid, near extremal, spin. These $M_\bullet$ and $a_\bullet$ values rule out both vector bosons and axions of masses $\sim10^{-16}$ eV. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16936v1

Evolutionary Origins of Binary Neutron Star Mergers: Effects of Common Envelope Efficiency and Metallicity by Monica Gallegos-Garcia et al. on Wednesday 30 November The formation histories of compact binary mergers, especially stellar-mass binary-black hole mergers, have recently come under increased scrutiny and revision. In this paper we revisit the question of the dominant formation channel and efficiency of forming binary neutron-star mergers. We use the stellar and binary evolution code MESA and implement an up-to-date and detailed method for common envelope and mass transfer. We preform simulations for donor masses between 8-20 solar masses with a neutron star companion of 1.4 and 2.0 solar masses, at two metallicities, using varying common envelope efficiencies, and two prescriptions for electron-capture supernovae. In contrast to the case of binary-black hole mergers, for a neutron star companion of 1.4 solar masses, all binary neutron star mergers are formed following a common envelope phase, while for a neutron star mass of 2.0 solar masses we identify a small subset of mergers following only stable mass transfer if the neutron star receives a large natal kick. Regardless of neutron star companion mass, we find that large supernova natal kicks are favored in the formation of binary neutron star mergers, and find more mergers at subsolar metallicity compared to solar. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15693v1

Anisotropic Photon and Electron Scattering without Ultrarelativistic Approximation by Anderson C. M. Lai et al. on Wednesday 30 November Interactions between photons and electrons are ubiquitous in astrophysics. Photons can be down scattered (Compton scattering) or up scattered (inverse Compton scattering) by moving electrons. Inverse Compton scattering, in particular, is an essential process for the production of astrophysical gamma rays. Computations of inverse Compton emission typically adopts an isotropic or an ultrarelativistic assumption to simplify the calculation, which makes them unable to broadcast the formula to the whole phase space of source particles. In view of this, we develop a numerical scheme to compute the interactions between anisotropic photons and electrons without taking ultrarelativistic approximations. Compared to the ultrarelativistic limit, our exact results show major deviations when target photons are down scattered or when they possess energy comparable to source electrons. We also consider two test cases of high-energy inverse Compton emission to validate our results in the ultrarelativistic limit. In general, our formalism can be applied to cases of anisotropic electron-photon scattering in various energy regimes, and for computing the polarizations of the scattered photons. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15691v1

A very luminous jet from the disruption of a star by a massive black hole by Igor Andreoni et al. on Wednesday 30 November Tidal disruption events (TDEs) are bursts of electromagnetic energy released when supermassive black holes (SMBHs) at the centers of galaxies violently disrupt a star that passes too close. TDEs provide a new window to study accretion onto SMBHs; in some rare cases, this accretion leads to launching of a relativistic jet, but the necessary conditions are not fully understood. The best studied jetted TDE to date is Swift J1644+57, which was discovered in gamma-rays, but was too obscured by dust to be seen at optical wavelengths. Here we report the optical discovery of AT2022cmc, a rapidly fading source at cosmological distance (redshift z=1.19325) whose unique lightcurve transitioned into a luminous plateau within days. Observations of a bright counterpart at other wavelengths, including X-rays, sub-millimeter, and radio, supports the interpretation of AT2022cmc as a jetted TDE containing a synchrotron "afterglow", likely launched by a SMBH with spin $a \gtrsim 0.3$. Using 4 years of Zwicky Transient Facility (ZTF) survey data, we calculate a rate of $0.02 ^{+ 0.04 }_{- 0.01 }$ Gpc$^{-3}$ yr$^{-1}$ for on-axis jetted TDEs based on the luminous, fast-fading red component, thus providing a measurement complementary to the rates derived from X-ray and radio observations. Correcting for the beaming angle effects, this rate confirms that about 1% of TDEs have relativistic jets. Optical surveys can use AT2022cmc as a prototype to unveil a population of jetted TDEs. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16530v1

Multiplicity of TeV muons in air showers detected with IceTop and IceCube by Stef Verpoest. on Wednesday 30 November The IceCube Neutrino Observatory at the South Pole can provide unique tests of muon production models in extensive air showers by measuring both the low-energy (GeV) and high-energy (TeV) muon components. We present here a measurement of the TeV muon content in near-vertical air showers detected with IceTop in coincidence with IceCube. The primary cosmic-ray energy is estimated from the dominant electromagnetic component of the air shower observed at the surface. The high-energy muon content of the shower is studied based on the energy losses measured in the deep detector. Using a neural network, the primary energy and the multiplicity of TeV muons are estimated on an event-by-event basis. The baseline analysis determines the average multiplicity as a function of the primary energy between 2.5 PeV and 250 PeV using the hadronic interaction model Sibyll 2.1. Results obtained using simulations based on the post-LHC models QGSJet-II.04 and EPOS-LHC are presented for primary energies up to 100 PeV. For all three hadronic interaction models, the measurements of the TeV muon content are consistent with the predictions assuming recent composition models. Comparing the results to measurements of GeV muons in air showers reveals a tension in the obtained composition interpretation based on the post-LHC models. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16970v1

Identifying the physical origin of gamma-ray bursts with supervised machine learning by Jia-Wei Luo et al. on Wednesday 30 November The empirical classification of gamma-ray bursts (GRBs) into long and short GRBs based on their durations is already firmly established. This empirical classification is generally linked to the physical classification of GRBs originating from compact binary mergers and GRBs originating from massive star collapses, or Type I and II GRBs, with the majority of short GRBs belonging to Type I and the majority of long GRBs belonging to Type II. However, there is a significant overlap in the duration distributions of long and short GRBs. Furthermore, some intermingled GRBs, i.e., short-duration Type II and long-duration Type I GRBs, have been reported. A multi-wavelength, multi-parameter classification scheme of GRBs is evidently needed. In this paper, we seek to build such a classification scheme with supervised machine learning methods, chiefly XGBoost. We utilize the GRB Big Table and Greiner's GRB catalog and divide the input features into three subgroups: prompt emission, afterglow, and host galaxy. We find that the prompt emission subgroup performs the best in distinguishing between Type I and II GRBs. We also find the most important distinguishing feature in prompt emission to be $T_{90}$, hardness ratio, and fluence. After building the machine learning model, we apply it to the currently unclassified GRBs to predict their probabilities of being either GRB class, and we assign the most probable class of each GRB to be its possible physical class. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16451v1

The Birth of a Relativistic Jet Following the Disruption of a Star by a Cosmological Black Hole by Dheeraj R. Pasham et al. on Wednesday 30 November A black hole can launch a powerful relativistic jet after it tidally disrupts a star. If this jet fortuitously aligns with our line of sight, the overall brightness is Doppler boosted by several orders of magnitude. Consequently, such on-axis relativistic tidal disruption events (TDEs) have the potential to unveil cosmological (redshift $z>$1) quiescent black holes and are ideal test beds to understand the radiative mechanisms operating in super-Eddington jets. Here, we present multi-wavelength (X-ray, UV, optical, and radio) observations of the optically discovered transient \target at $z=1.193$. Its unusual X-ray properties, including a peak observed luminosity of $\gtrsim$10$^{48}$ erg s$^{-1}$, systematic variability on timescales as short as 1000 seconds, and overall duration lasting more than 30 days in the rest-frame are traits associated with relativistic TDEs. The X-ray to radio spectral energy distributions spanning 5-50 days after discovery can be explained as synchrotron emission from a relativistic jet (radio), synchrotron self-Compton (X-rays), and thermal emission similar to that seen in low-redshift TDEs (UV/optical). Our modeling implies a beamed, highly relativistic jet akin to blazars but requires extreme matter-domination, i.e, high ratio of electron-to-magnetic field energy densities in the jet, and challenges our theoretical understanding of jets. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16537v1

Inferencing Progenitor and Explosion Properties of Evolving Core-collapse Supernovae from Zwicky Transient Facility Light Curves by Bhagya M. Subrayan et al. on Wednesday 30 November We analyze a sample of 45 Type II supernovae from the Zwicky Transient Facility (ZTF) public survey using a grid of hydrodynamical models in order to assess whether theoretically-driven forecasts can intelligently guide follow up observations supporting all-sky survey alert streams. We estimate several progenitor properties and explosion physics parameters including zero-age-main-sequence (ZAMS) mass, mass-loss rate, kinetic energy, 56Ni mass synthesized, host extinction, and the time of explosion. Using complete light curves we obtain confident characterizations for 34 events in our sample, with the inferences of the remaining 11 events limited either by poorly constraining data or the boundaries of our model grid. We also simulate real-time characterization of alert stream data by comparing our model grid to various stages of incomplete light curves (t less than 25 days, t less than 50 days, all data), and find that some parameters are more reliable indicators of true values at early epochs than others. Specifically, ZAMS mass, time of explosion, steepness parameter beta, and host extinction are reasonably constrained with incomplete light curve data, whereas mass-loss rate, kinetic energy and 56Ni mass estimates generally require complete light curves spanning greater than 100 days. We conclude that real-time modeling of transients, supported by multi-band synthetic light curves tailored to survey passbands, can be used as a powerful tool to identify critical epochs of follow up observations. Our findings are relevant to identify, prioritize, and coordinate efficient follow up of transients discovered by Vera C. Rubin Observatory. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15702v1

Studying cosmic ray sources using intergalactic electromagnetic cascades by Anna Uryson. on Wednesday 30 November In this paper intergalatic electromagnetic cascades are used as a probe of cosmic ray sources. This is achieved as follows. In extragalactic space cosmic rays initiate electromagnetic cascades in which gamma-ray and neutrino emission arises. We used the joint analysis of cosmic ray data, along with extragalactic gamma-ray and neutrino emission, to study particle acceleration in the vicinity of supermassive black holes. Particle injection spectra depend on processes of particle acceleration, and here we discuss models with various injection spectra. The computation of the propagation of cosmic rays in space were performed using the publicly avaliable TransportCR code. It was found that a new subclass of sources might exist that does not contribute to the particle flux on Earth, instead to gamma-ray and neutrino emissions arising in electromagnetic cascades. In addition, the upper limit of the relative number of 'exotic' supermassive black holes surrounded by a superstrong magnetic field is derived. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.16556v1

Observing Signals of Spectral Features in the Cosmic-Ray Positrons and Electrons from Milky Way Pulsars by Ilias Cholis et al. on Wednesday 30 November The Alpha Magnetic Spectrometer (AMS-02) has provided unprecedented precision measurements of the electron and positron cosmic-ray fluxes and the positron fraction spectrum. At the higher energies, sources as energetic local pulsars, may contribute to both cosmic-ray species. The discreteness of the source population, can result in features both on the positron fraction measurement and in the respective electron and positron spectra. For the latter, those would coincide in energy and would contrast predictions of smooth spectra as from particle dark matter. In this work, using a library of pulsar population models for the local part of the Milky Way, we perform a power-spectrum analysis on the cosmic-ray positron fraction. We also develop a technique to cross-correlate the electron and positron fluxes. We show that both such analyses, can be used to search statistically for the presence of spectral wiggles in the cosmic-ray data. For a significant fraction of our pulsar simulations, those techniques are already sensitive enough to give a signal for the presence of those features above the regular noise, with forthcoming observations making them even more sensitive. Finally, by cross-correlating the AMS-02 electron and positron spectra, we find an intriguing first hint for a positive correlation between them, of the kind expected by a population of local pulsars. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15709v1

Fundamental physics with neutron stars by Joonas Nättilä et al. on Wednesday 30 November Neutron stars are rich laboratories of multiple branches of modern physics. These include gravitational physics, nuclear and particle physics, (quantum) electrodynamics, and plasma astrophysics. In this chapter, we present the pioneering theoretical studies and the pivotal historical observations on which our understanding of neutron stars is based on. Then, we discuss the usage of neutron stars as probes of fundamental theories of physics. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15721v1

Hard X-ray Observations of the Hydrogen-poor Superluminous Supernova SN 2018hti with NuSTAR by Igor Andreoni et al. on Wednesday 30 November Some Hydrogen-poor superluminous supernovae are likely powered by a magnetar central engine, making their luminosity larger than common supernovae. Although a significant amount of X-ray flux is expected from the spin down of the magnetar, direct observational evidence is still to be found, giving rise to the "missing energy" problem. Here we present NuSTAR observations of nearby SN 2018hti 2.4y (rest frame) after its optical peak. We expect that, by this time, the ejecta have become optically thin for photons more energetic than about 15keV. No flux is detected at the position of the supernova down to $F_{\rm{10-30keV}} = 9.0\times 10^{-14}$ erg cm$^{-2}$ s$^{-1}$, or an upper limit of $7.9 \times 10^{41}$ erg s$^{-1}$ at a distance of 271Mpc. This constrains the fraction of bolometric luminosity from the putative spinning down magnetar to be $f_{\rm X} \lesssim 36$% in the 10-30keV range in a conservative case, $f_{\rm X} \lesssim 11$% in an optimistic case. arXiv: http://arxiv.org/abs/http://arxiv.org/abs/2211.15749v1